function [S] = plotspeed(stop, axe, particles)
    S = zeros(size(particles,2),stop);
    for i = 1:1:stop
        [name,err] = sprintf('speed_%d.dat', i);
        A = load(name);
        S(1:1:size(particles,2), i) = A(particles, axe);
    end
